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Abstract: The aim of this paper is to estabhsh a relationship between the number of 
galaxies and the rate of type la supernovae. In addition, present approach prepares a 
framework to analyze this relationship and its dependence on the redshift (z), giving aid 
to the deepening of studies on radial distribution of galaxies. The work also shows how 
the counts interrelates to topological aspects regarding the evolution of the universe from 
z ^ 2. The relationship was examined apart from morphological considerations. 
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Behold the heavenly space, its stability, celerity, and 
quit, finally, to look vile things. 

Boecio 

The existing type la supernovae scenario 

The supernovae counting 

Type la supernovae are among the biggest 
thermonuclear explosions in the Universe. During the 



two last decades, they gained great importance in 
cosmology as benchmarks to measure distances due to 
their regular behavior like standard candles and because 
of their essential role in the search for decisive facts to 
corroborate the accelerated expansion of the Universe. 
At first, there was a belief that type la supernovae 
(SNe) originate only from a carbon-oxygen white dwarf 
(WD) in a binary system, where the dwarf companion 
is a red giant. Nowadays, two possible progenitor 
models are discussed, that is, a) one starting from the 
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classical hypothesis that the WD, accreting from its 
main sequence or red giant companion, grows in mass 
until it reaches the critical Chandrasekhar limit, and 
b) other starting from the hypothesis of a close double 
WD system merging after orbital shrinking due to the 
emission of gravitational wave radiation. In both cases, 
the nature and evolution of the binary system remains 
a puzzle. The time elapsed from the birth of the binary 
system to the SN explosion, called "delay time", spans 
from tens of millions of years to ten billion years or even 
more. As a consequence, the counting of SNe la (or the 
rate of generation) reveals the star formation history of 
a galaxy accordingly the distribution of the delay times. 
The analysis of the SNe la rate as a function of redshift 
and the inquiry whether their properties evolve with 
redshift would be powerful tools to aid the investigation 
of the nature of the galaxy environments and their 
connections with the rate of galaxy formation itself. 
Type la supernovae have been key pieces in studies 
on the true nature of the distribution of matter in the 
universe. The large number of uncertainties referred to 
the immensity of the distances and time scales leaded 
to alternatives to the more accepted interpretation 
of the data, which was taken in accordance with the 
homogeneous Friedmann-Lemaitre- Robertson- Walker 
(FLRW) cosmology. In this way, a simple and realistic 
option is provided by the inhomogeneous Lemaitre- 
Tolman (LT) cosmology, applied and studied by several 
authors. For instance, there is an article by D. Garfinkle 
in which an inhomogeneous LT spacetime is treated 
in an understandable and very p ragmatic way t o 
modeling type la supernovae data (iGarfinkld . |2006[ ). 
By its simplicity and clarity, I consider this study as 
a fundamental starting point for further research on 
inhomogeneous cosmology. Also based on the adoption 
of a LT cosmology, K. Bolejko analyzes whether the 
modeling of inhomogeneous distribution of matter with 
no cosmological constant is realistic or not, taking as 
object of investigati on the r educt ion of brightness of 
type la supernovae (|Boleikd . 120081 ). The SNe rate at a 
given redshift is usually computed as the ratio between 
the number of known SNe la and the control time of 
the monitored galaxies at that redshift. One may object 
that statistics would be a poor tool to model a tracer 
for the rates both at low redshift, due to the difficulty 
of sampling large volumes, and at high redshift, due to 
the difficulty of detecting and typing faint SNe. However 
they are still seen as rare and transient events, thousands 
of SNe la are now discovered per year (jDecadal, 2010), 
and the deep sky observations with suitable time interval 
required to enlarge the number of real SNe candidates 
are more and more attainable as technology advances. 
Thus, the more the range of z is wide, the more will 
be the reliability on the proposed z-dependence of the 
rate. The main motivations of this work were a)- to 
relate quantitatively the countings of galaxies and type 
la supernovae, trying to bring new elements to contrast 
conventional searching on radial distribution of galaxies, 
b)- to find and discuss possible evolutionary topological 



properties of the Universe from z ^ 2 based on that 
quantitative relation, and c)- to question weather the 
eventual existence of evolutionary "jumps" in the 
Universe would provoke relevant anisotropics on CMB. 
The central idea is to define a relation between the 
number of galaxies and the SNe rate by redshift, since 
we know that the energy deposition from supernovae 
into the environment is substantial to star formation and 
even to galaxy formation. In addition, discounting all 
troubled observational biases, such as that caused by the 
interstellar dust, which produces more extinction of the 
high redshift supernovae, distinct rates for different high 
redshifts could result from the evolution of progenitor 
systems by variation in the mass, composition, and 
metallicity of type la supernovae. The test of that 
hypothesis requires accurate measurements of SN rates 
at various cosmic epochs. Beside, to minimize the 
uncertainty in the estimates of SN rates, a statistically 
significant SNe sample is necessary. The chemical 
enrichment of the galaxies by the returns of synthesized 
elements from the explosive burning processes of 
supernovae is on its own a strong reason to investigate 
the explosion rates by redshifts and this is a matter of 
many works. The literature on galaxy and supernova 
countings is copious. At first, I distinguish the works 
on supernova rate by Pain et al (1996), Cappellaro et 
al (1999) and from the EROS CoUaboration (2000), 
followed by Blanc et al (2004) and culminating with 
Sharon et al (2010). In particular, about galaxy number 
counts and clusters of galaxies we have fine works from 
Campos (1995), Metcalfe et al (2000), Grahan et al 
(2008) and Labini et al (2009). Another fundamental 
study is the classical paper of Lilly et al (1995). As low 
redshifts are the regions where the effective work to 
understanding SNe la must be done, Cappellaro et al 
proposed to estimate the SNe rate in the local Universe 
as a function of both galaxy morphological type and 
colors. On the other hand, the EROS Collaboration, 
searching nearby supernovae from z 0.02 to 0.2, took 
the explosion rate of J^obs observed type la supernovae 
by the expression 



n = Afobs/S, 
where 



gal i 



(1) 



+ 00 



: {t,Zi)dt, 



(2) 



with {t, Zi) as the efficiency to detect in galaxy i a 
type la supernova whose maximum occurs at time t in 
the supernova rest frame. Galaxy i has a blue luminosity 
Li, and the integral {t, Zi) dt is called its "control 

time" . For the EROS simulations, the sum S was proved 
to be equal to 9.09 x 1O^^/i~^L0b yr and the error is 
statistical and usually given at a 68% confidence level. 
Further, also according to EROS collaboration, there 
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must be a cut on the galaxy magnitude, so that we 
may ensure that the host galaxy for a type la supernova 
is classified as such during the visual scanning. For 
example, at z = 0.1, the cut corresponds to an absolute 

magnitude of Mgai 20. + 5log{Ho/60kms~^Mpc'^) 

for the host galaxy. Possible samples of supernovae are 
available in some good catalogs. An interesting sample 
from ESSENCE Supernova Survey (Miknaits+, 2007) is 
shown in Table 1 with probabilities for the real type la 
supernovae candidates. In this paper, I present a general 
form to predicting the number of galaxies associated to 
the type la supernova explosion rate. I assumed that 
the expected rate at a given red shift is proportional to 
the amount of galaxies at the redshift. In this approach 
there are no concerns about morphology but purely 
statistical considerations. Since in this study no attempt 
was made to derive rates for different galaxy types and 
the individual SN rate was implicitly assumed depending 
on host galaxy luminosity, the predicted number of 
galaxies is such that the rate is achieved by the presence 
of enough galaxies with appropriated luminosity. 
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Frameworks to counting galaxies from SNe 
explosion rate 

The classical galaxy counting 

The numerical differential counting of objects (galaxies 
or clusters, depending on the numerical density and the 
source masses), or solely differential counting, is defined 
in terms of the volumn as 



dN = n{t)dV, 



(3) 



where n{t) is the numerical density {V ^) and dV the 
volumn element. This last is givem by 



dV = ^di 



: sin OdOdcf), 



dV 



Aiva^r'^dr 
Vl - kr'^ ' 



(3.a) 



(3.&) 



where a is the scale factor or expansion rate of the 
Universe, r is the commoving radius. A; is a constant of 



curvature and g is the determinant of the Robertson- 
Walker metric tensor. From these expressions we deduce 
that dV is a spherical shell with radius ar and thickness 
(1 — kr'^)^^/^adr. In terms of redshift, the connection 
between the forecasted behavior of the differential 
counting and that builded from the observations is given 

by 



dN 
dz 



obs 



ij{z) dN 
nc{z) dz ' 



(4) 



where ipi^) is the selection function applied on the 
required catalog and nc{z) is the radial numerical density 
in terms of the commoving volumn. A more complete 
approach on this subject may be found, for instance, in 
Iribarrem (2009). 

First approach 

Let us take the EROS surveying; beside, let us assume 
that the star formation activity is s t eady and continuous 



for most of galaxies (jShimasakul . I1997I ). and for any 
redshift, starting from the requirement that one SN must 
be located at one host galaxy, we have a SN rate and a 
differential galaxy counting related by 



n = T(z) 



dN 
dz 



S 



(5) 



which means there is a relation between the number of 
galaxies and the number of type la supernovae, strongly 
depending on z, with T(z) being a certain function of the 
redshift to be defined ahead. Since dN/dz is independent 
of the cosmology, getting the differential counting for a 
certain catalog we may choose the cosmological model 
further. From equations (1) and (2) we have 



dNobs = A/'ofc. 



dz 



T(z)5' 



Nobs = / J^ob. 



dz 



'T(z)5' 



Nobsdz 



T(^) E UC^^^{t,Z.,)dt 
gal i 



(6) 
(7) 

(8) 



As Zi refers to the i redshifts of a certain observed galaxy 
sample, the differential dz takes account only of the 
integration of 1/T(z). For S = 9.09 x IQ^'^h-'^LQB yr, 
equation (8) reduces to 



Nobs = 



Nobsdz 



T{z) X 9.09 X W^Ii-'^Lqb yr 



(9) 



According to EROS, there were surveyed 80 square 
degrees, comprising eight detected supernovae, four of 
which spectroscopically identified as type la supernovae. 
Thus, 



dz 



9.09 X 1012/1-2^05 yr J T(z) 



(10) 
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The general host function 

On the other hand, accordingly Pain et al (1996) in an 
interesting study on the expected number of supernovae 
as a function of z, the observed SN rate is also given by 



7^ = 



(1 + ^)' 



(11) 



where r^n is the rate in the rest-frame of the supernovae. 
Thus, we have 



T(z) 



dN 



dz 



Mob 



obs 
Nobs 



s ' 



(1 + ^) S ' 
From equations (12) and (13) we get 



dN 



dz 



= 1/T(z)-^, 

obs (1 + 



T(z)(l + z)- 
Now we equal equations (8) and (15), so that 



Nobsdz 



T(^) E Uj:^e,{t,z,)dt 

gal i 

J^obs 



idz 



T(z)(l + z) 



nz)T. L,^^^^{t.^^)dt T(Z)(1 + Z)' 
gal i 



Mob 



(1+^) 



(t, Zi) dt, 



(12) 
(13) 

(14) 
(15) 

, (16) 
(17) 

(18) 



which is precisely equation (13). The function ^ = 
fsji/T(2;)(l + z) is called "host function" and it brings 
the cosmological phase signature at a given z for a certain 
host galaxy. The signatures are positive till z ^ 0.5 
(figure 1). The change in the signal of the host function is 
consistent with the change from a decelerating Universe 
(matter dominated phase) up until about 5 billion years 
ago (around z = 0.5) to the current accelerating stage 
( Lidmanl - liool . At certain z, the integration of the host 
function gives the expected number of galaxies related to 
the expected rate of supernovae; so, function T(z) acts 
upon the derivative of the number of galaxies at that z to 
give the rate of supernovae for the respective z-shell. As 
very rare events, statistically speaking, we may suppose 
from the viewpoint of the observer one supernova to one 
galaxy, so that the rate of supernovae is in sum closely 
the rate of host galaxies. However, as pointed out by 
Broadhurst, Taylor and Peacock (1994), we may assume 
some cosmological premises to fit the data of redshift 
surveys, remembering that the representative function 
of the number of galaxies Nobs is model- independent. 
The numerical results for Nobs were obtained from the 
given luminosity function and these authors worked at 
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Figure 1; evolution of host function ^ (grey) and number of 
galaxies Nobs (blue) until 2 = 2. 
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Figure 2: evolution o/T (blue) and ^ (grey) until z = 2. 

the faint limit for spectroscopy in the R band, in which 
R — 22.5. They adopted for this case the expression 



Nobs = 11.7z^-''^ea;p[-(z/0.51) 



1.791 



(19) 



where Nobs refers to the probability distribution for 
redshift; I used that to infering the host function 
precisely for this case. Equations (15) and (16) furnish 



11.7zi-«3exph(./0.51)i-^^] ^ I ^ 
from which we gain 



dz 



(z)(l + z)' 



(20) 



T(z)(l + z) 



19.07zO-63(l - 3.66zi-79)g-.i ™/o.2996 (21) 
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The last equation gives 



VO.2996 



19.07^0-63(1 -3.66^i-79)(l + z) ■ 



(22) 



Substituting this result, the integration in equation (15) 
furnishes 



N, 



obs 
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detect some additional harsh distortion in the expected 
gained energy of the electrons that, crossing the core 
of a galaxy cluster, are Thomson scattered distorting 
the CMB blackbody spectrum. The essential point to be 
investigated is the way one might decide on the origins of 
certain anisotropy in the cosmic microwave background. 
For instance, the spectral shapes of the anisotropy 
generated by the gravitational lensing are identical to 
those generated by kinematic Sunyaev-Zel'dovich effect. 



(-50 2^+183^^)^^ + 1. (23) 

We note that T has a discontinuous structure with 
different limits at a particular point around 2; ~ 0.5 
depending on whether this point is approached from the 
positive or negative directions (Figure 2). In addition, 
the two "foils" of T are not symmetric. In fact, one may 
define a map Ai / : 1(2) =^ ^ to cover T, that is, to 
bypass the discontinuity. This covering sounds artificial, 
and I defend the opinion that the T{z)-iump a,t z ^ 0.5 
must be a real break in the cosmic evolution related 
to the harsh phase transition toward an accelerating 
Universe (the hardness to understand or to explain the 
physics of the jump is not acceptable as an argument 
to bypass discontinuity). The discontinuous structure 
T is featured by one singularity. So, since limits are 
topological constructions, that is, manipulations with 
limits are based on the topology of number spaces, let z 
be a singular point around z 0.5; also let := {Z,T) 
be a topological space. We say that z is an isolated 
point of 3 := (Z, T) ^-.^^ {z} G T, such that all other 
points are limit points enclosed in the set Lim{^) := 
{z Cz Z \ z ^ T}. In this case, 3 is said to be a singular 
space. Thus, defining the function T{z) we are building 
a map : =^ 3' that carries topological space 
into topological space 9'. Both spaces are subspaces 
of M. The intervals of z arc < z < ¥ and z < z < 
2; the union of the two non-symmetric disjoint foils 
gives T. Because of its singularity, which means that 
we need to manipulate limits, and by the fact that it 
is not depending on the metric (cosmology), function 
T{z) was called "Topological Evolution from High 
Redshift Astronomy" (TEHRA). This function may 
bring some interesting additional insights to relativistic 
cosmology; as partial differential equations, Einstein's 
equations describe only local properties of spacetime 
and do not give any information about the global 
structure (the topology) of the spacetime. In physics, 
topological effects connected to non-genetic holes or 
harsh breakouts on manifolds are generally expensive in 
energy. This cost in energy must be detectable such as 
anisotropics on CMB. Nevertheless, in present state of 
knowledge it is very difficult to decide if the TEHRA 
singularity (TS) originated some kind of anisotropy 
on CMB. Such anisotropy, if detected, would be an 
indirect evidence of TS after an accurate analysis of 
all possibilities. Recently, I suggested to my colleagues 
a project to investigate exhaustively the Sunyaev- 
Zel'dovich effect from z < 0.5 and z > 0.5 trying to 



Conclusion 

Present work explained what I called "host function" as 

an attempt to connect type la supernovae rate and the 
counting of galaxies in accordance to the observational 
data. The central idea is that for a certain rate of 
supernovae there are sufficient host galaxies to warrant 
the rate. There is a wide field for detailed research in 
this way to discussing the host function. For instance, 
some possible physical consequences of evolutionary 
topological features pointed by the model were 
explained, suggesting further investigations and tests on 
CMB to look for relics from the singularity predicted. In 
addition, however the counting of galaxies is independent 
on cosmology, it would be particularly interesting to take 
into account the inhomogeneous Lemaitre-Tolman and 
the Friedmann-Lemaitre-Robertson- Walker cosmologies 
to compare the distribution of galaxies deduced from the 
above theoretical approach with that predicted by both 
cosmological models. 



References 

Blanc, G. et al. (2004), Type la supernova rate at a 
redshift of 0.1 , arXiv:astro-ph/0405211. 

Bolejko, K. (2008), Supernova la observations in the 
Lemaitre-Tolman model, arXiv:astro-ph/0512103. 

Broadhurst, T., Taylor, N. & Peacock, J. (1994), 
Mapping cluster mass distributions via gravitational 
lensing of background galaxies, arXiv:£istro- 
ph/9406052. 

Campos, A. (1995), Galaxy number counts, astro- 
ph/9510051. 

Cappellaro, E. et al. (1999), A new determination of 
supernova rates and a comparison with indicators for 
galactic star formation, arXiv:astro-ph/9904225. 



6 Nilo Serpa 



Decadal Survey of Astronomy and Astrophysics (2010- 
2020), Type la supernova science 2010-2020, white 
paper submitted to the Astro2010. 

Garfinklc, D. (2006), Inhomogcncous spacctimcs as a 
dark energy model, arXiv: gr-qc/0605088. 

Grahan, M. et al. (2008), Type la supernovae rates and 
galaxy clustering from the CFHT Supernova Legacy 
Survey, arXiv:0801.4968. 

Iribarrem, A. (2009), Contagem relativfstica de fontes 
cosmologicas e a fungao de luminosidade galactica, 
Magister Thesis, UFRJ - BR. 

Labini, F. et al. (2009), Large-scale fluctuations in the 
distribution of galaxies from the Two Degree Field 
Galaxy Redshift Survey, arXiv:0902.0229. 

Lidman, C. (2004), Observing distant type la supernovae 
with the ESQ VLT, ESO - Reports from 
Observers. 

Lilly, S. et al. (1995), The evolution of field galaxies. 
International Astronomical Union, Apj, 455, 
108. 

Metcalfe, N. et al. (2000), Galaxy number counts - 
V. ultra-deep counts: The Herschel and Hubble deep 
fields, arXiv:astro-ph/0010153. 

Perlmutter, S. et al. (1999), Constraining dark energy 
with type la supernovae and large-scale structure, 
Astrophys. J. 517 565. 

Pain, R. et al. (1996), The type la supernova rate at 
z ^ 0.4, The Supernova Cosmology Project: IV. 

Sharon, K. et al. (2010), The type la supernova rate in 
redshift 0.50.9 galaxy clusters, arXiv:1006.3757. 

Shimasaku, K. & Fukugita, M. (1997), The history of 
galaxies and galaxy number counts, arXiv:astro- 
ph/9705057. 

The EROS collaboration. Type la supernova rate at z ~ 
0.1, Astron. Astrophys. 362, 419-425. 



